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ABSTRACT 


The goal of this thesis is to evaluate the use of surface waves to detect buried objects. The 
source used to generate the surface waves was a three element phased array controlled by 
Lab VIEW visual instruments. This research included developing a source and receiver, evaluating 
attenuation and azimuthal dependence of the surface wave propagation, detection using scattering 
within a tank of sand, and using the three element array to beamform surface waves. It was 
successfully demonstrated that target localization using surface wave scattering and beamforming 
with a phased array is possible. 
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I. INTRODUCTION 


A. OBJECTIVE 

The purpose of this research was to evaluate the feasibility of 
using a phased array to generate surface waves to detect buried 
land/sea mines. 

B. MOTIVATION 

The detection of buried mines in the ocean bottom, surf zone 
and on land is a continuing civil and military problem. In October, 
1994 the United Nations Secretary General Boutros Boutros-Ghali 
commented that the estimated number of mines still in place from 
conflicts of this century is over 100 million through within 62 
countries [Ref 1]. These mines, originally intended to disrupt 
commercial activity, utility plant operations, port facilities, 
transportation routes, and troop movement take their toll on 
postwar civilian populations throughout the world. The technology 
to safely and efficiently remove land mines has not been developed. 
The current cost of removing a single $25 anti-personnel mine is 
estimated by the United Nations [Ref 1] to ranges from $300 - 
$1000. This cost is projected to be higher as popular plastic mines 
render magnetic mine clearance devices currently in use obsolete. 

The strategic advantage gained from the Iraqi mine field in the 
Persian Gulf can be debated, given the results of the war. However, 
there is no question that it was effective in damaging two US. 
warships, altering the assault objective to Faylakah Island, and 
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preventing the amphibious assault on Kuwait. At the end of the 
Persian Gulf war 1,300 mines were recovered with 24% being 
bottom mines in various states of scouring according to the Center 
for Security Strategies and Operations [Ref 2]. These mines could 
have had a much greater impact on the outcome of the war and the 
number of personnel and material losses had a different assault plan 
been chosen. Set up as a barrier to amphibious assault, the Iraqi 
mine field was successfully avoided and removed once it was 
discovered. If the mines had been used more effectively to hamper 
the positioning of ships or if amphibious assault had been an only 
option, mine clearance would have been a much bigger problem than 
as it turned out. 

As implied in “. . .From the Sea,” for our current naval 
doctrine to work, some method of rapid detection of mines must be 
developed to deal with buried mines. 

C. BACKGROUND 

For more than 25 years the US. Navy has utilized marine 
mammals to detect objects in the water and buried in the ocean 
floor. The continuing study of these animals and how their sensor 
systems work could provide valuable insights for the next 
generation of sonar for mine detection. In a similar fashion, sand 
scorpions of the Mojave desert, who are able to capture their prey by 
detecting Rayleigh waves, might have some of the answers for 
buried mine detection, [Ref 3]. Utilizing their eight independent 
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tarsal slit sensilla as vibration receptors, scorpions are able to 
detect small vibrations in the sand resulting from the Rayleigh 
surface waves and compression waves (P waves) created by the 
movement of their prey. They are then able to accurately estimate 
time differences between arrivals of the two types of waves in 
order to formulate target distance and direction from the timing and 
ordering of signals received by their sensors. At distances of 10 cm 
they can derive target angle and distance very accurately and tests 
have shown they can detect disturbances out to 30 cm [Ref 3]. 

In elastic media, Rayleigh waves propagate with a non- 
dispersive speed slower than the shear and compression wave 
speeds of the medium. Because Rayleigh waves are localized to 
within a few wavelengths of the surface of the medium they could 
be ideal for hunting mines that must be placed near the surface to be 
effective. From 1987 to 1992, BBN Systems and Technology 
undertook a feasibility study to design an experimental system for 
detecting and locating buried land mines using acoustical probing 
[Ref 4]. In their report, BBN employed source and receiver arrays 
which achieved detection distances of 15 feet. The source used in 
that work was a water filled tube in which a high voltage electrode 
created a spark. The resulting pressure wave was transmitted 
through a rubber membrane in contact with the ground. It is 
estimated that only 67% of the energy from a point source energy is 
propagated in the form of Rayleigh waves. The remaining energy is 
disbursed in the form of downgoing vertically polarized shear (SV) 
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waves (23%) and compression (P) waves (10%). Coupled with no 
source directivity, BBN's system is inefficient as a method of 
generating surface waves. 

The focus of the research presented in this thesis is to build 
an improved surface wave source. The source will consist of a 
phased array of several discrete sources, operated at an optimal 
frequency for the specific soil type. 

D. STEPS OF RESEARCH 

The stages of research involved: 

1. Building an appropriate test tank. 

2. Developing sources and receivers. 

3. Develop data acquisition software using LabVIEW. 

4. Measuring the performance of a three element array of 
sources. 
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II. DEVELOPMENT AND DESIGN 


The purpose of this chapter is to 1) review the characteristics 
of Rayleigh waves, 2) discuss the test tank, source and receiver used 
in this research, and 3) discuss the LabVIEW-based data acquisition 
system. 

A. WAVE PROPAGATION IN SOLIDS 

A point source acting at the surface of an elastic medium 
excites three basic types of waves: compressional waves, shear 
waves, and surface waves. Compressional waves are characterized 
by particle motion along the direction of propagation, whereas in 
shear waves the motion is transverse to the direction of 
propagation. Both of these types of waves can exist within the bulk 
of the solid and spread spherically away from a point source. 

Surface waves, as the name implies, are confined to the surface of 
the medium [see Figure (1.1)]. Their amplitude decays exponentially 
with distance from the surface, with a decay constant 
approximately equal to the wavelength. As a result of their 
confinement to the surface, the propagation is cylindrical in nature. 
The particle motion is a retrograde ellipse, confined to the vertical 
plane. Figure (1.2). The speed of propagation is less than the shear 
wave speed for the medium. Surface waves are commonly referred 
to as Rayleigh waves, although strickly speaking, Rayleigh wave is 
the correct name only for the case of a wave on the surface bounding 
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Figure (1.1) Compressional wave and Rayleigh wave. 



Figure (1.2). Rayleigh wave propagaton. 


a semi-infinite elastic medium and vacuum. However, for practical 
purposes the physics of Rayleigh waves accurately describes the 
practical situation of a wave on the surface of a solid, many 
surface-wavelengths thick, with air above. 










B. TEST TANK 

In this thesis we want to investigate the propagation of 
surface waves at an air/wet sand interface in a controlled, 
laboratory environment. The first step is to establish a suitable 
environment. Preliminary tests were performed in a 1.1 m x 3 m, 
fiberglass lined wooden tank, filled with approximately 60 cm of 
wet sand. The purpose of these preliminary measurements were to 
determine approximate surface wave speeds, and therefore 
wavelengths. The results indicated wave speeds in the range 60 - 
90 m/s. The corresponding wavelengths for frequencies from 300 - 
2000 Hz are approximately 25-4 cm. Therefore, we need a sand 
depth of approximately 1 m to mimic a semi-infinite medium at the 
lowest frequencies of interest. The width of the tank needed to be 
as large as possible to prevent interference with reflections from 
its sides. These conditions were approximately met by a twelve 
foot diameter, four foot deep redwood tank. This tank was filled 
with approximately 17 tons of medium gain (.6 mm) beach sand from 
Moss Landing in Monterey County, California. The sand was washed 
to remove salt and filtered for large objects. Using a garden hose 
with spray attachment the surface was uniformly wetted prior to 
each experiment. 

C. THE SOURCE 

The source design was guided by the following requirements: 

1. It should be a vertical point source that could generate 
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signals with sufficient amplitude to allow detection of 
surface waves anywhere in the test tank. 

2. It should provide a repeatable waveform. 

3. It should be capable of high repetition rates. 

4. It should have good ground coupling that would last for an 
extended period. 

5. It should provide clean signals over a range of frequencies, 

from approximately 300-2000 Hz. 

In addition, it had to be relatively inexpensive and easily 
modified given the time and resources available. It was therefore 
not the intent of this research to develop a working unit for mass 
production. After various types of drivers were tested it was found 
that a four inch diameter, high cone excersion, electrodynamic 
speaker with an aluminum 45 degree angle cone attached to the 
center of the diaphragm, worked well. In addition, a 25 gram mass 
was added to the cone appendage to improve ground coupling. 
(Figure(1.3), factsheet Appendix (A). 



Figure (1.3) Source design. 
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D. THE RECEIVER 


The design requirements for a receiver to detect ground 
vibration for this research were to provide: 

1. Adequate sensitivity to ground motion generated by the 

source within the distances of the tank. 

2. Broadband response. 

3. Good ground coupling that could give repeatable results. 

Several different designs were tried. The initial design 
focused on using a piezo-electric disk. The disks were acquired 
from small audio speakers available from any commerical 
electronics shop. The disks had appendages glued to them to provide 
coupling to the ground. They also had a mass attached with the top 
to improve coupling and reduce ringing. Figure (1.4) shows the 
different appendages that were tried. 


golf tee 


V 


flat based 
cylinder 



disk 
resting 
on surface 


Figure (1.4) Receiver ground coupling methods 
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The golf tee allowed for the greatest sensitivity however the 
received amplitude varied by greater than five percent. The inverted 
golf tee resting on the surface with approximately 20 grams of mass 
loading on a ring element reduced variations to within five percent. 
This coupled with signal averaging discussed later gave marginally 
repeatable results. However, the results were not considered 
satisfactory. Further tests were done using a commercial single 
axis geophone made by Halliburton Geophysical Servicves, model 
number S4408A. Though these geophones were not as sensitive as 
the various piezo-electric sensors, they were found to be more 
reliable in producing repeatable data. 

E. SUPPORTING EQUIPMENT 

The layout of the test tank and supporting equipment is shown 
in Figure (1.5). The setup consisted of five elements: 

1. The redwood tank filled with sand. 

2. Support equipment above the tank. 

3. Lab instruments including: 

a. An HP 3314A Function Generator which was used as a 
source function generator for attenuation measurements and to 
characterize cylindrical spreading. It was also used to provide a 
reference signal when designing the array function generator visual 
instrument (vi) created in LabVIEW. 

b. An COS 6100A Oscilliscope to visually monitor analog 
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input and output signals. 

c. An HP 3562A Dynamic Signal Analizer was utilized to 
isolate optimal propagation frequencies for given conditions. 

d. An Ithaco 1201 Low Noise Preamplifier was used to 
filter and amplify to the received signal. 

e. A Techron 5507 Power Supply amplified the source 


signal. 



4. A Macintosh llci computer with a LabVIEW interface board 

provided data aquisition. 

5. Sources and receivers 







F. DATA AQUISITION SYSTEM 

It was decided to perform data acquisition functions with 
LabVIEW. The data acquisition functions are implemented through 
software visual instruments or vis. The various vis used in this 
research are explained later and in the appendixes. The vis allow 
pulse generation, signal averaging, filtering, A/D conversion, and a 
high data recording rate. In addition to the software, two plug-in 
board were needed, the NB-MIO-16XL Data Aquisition board and the 
NB-AO-6 Analog Output board. 
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III. CHARACTERISTICS OF A POINT SOURCE 


The first step toward building a phased array surface wave 
source is to have well characterized single sources. The purpose of 
this chapter is to discuss the characterisitics of the waves 
generated by a single source. First, we examine the radial and 
azimuthal dependence the wave for different frequencies. Next, 
results are presented of scattering measurements from a buried 
target. Finally, we discuss target localization using multiple 
receiver, bistatic scattering. 

A. RADIAL DEPENDENCE 

The purpose of measuring the radial dependence of the wave 
propagation is threefold: 1) to determine if the medium is isotropic, 
2) to determine if the waves spread cylindrically, and 3) to 
determine the attenuation coefficient for later use in calculating 
beam patterns for the phased array. 

1. Procedure 

The point source was mounted on the beam in the middle of the 
tank after the sand had been wetted uniformly. Test points were 
laid out every 10 cm along two lines, oriented at right angles. The 
HP 3314A Function Generator was used to generate a single cycle 
pulse. This pulse was amplified with the Techron 5507 Power 
Supply and sent to the driver. The received signal was sent to the 
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the Ithaco 1201 Low Noise Preamplifier to be amplified and filtered. 
The signal was then recorded using the Intensity Field vi as 



shown in Figure (3.3). Measurements were taken from the outer 
radial distance moving toward the source in order not to disturb the 
sand associated with successive measurements. As can be seen 
from the graph, Figure (3.1), the signal decreases monotonically 
with range. To determine whether the spreading is cylindrical and 
the value at the attenuation coefficient (a), the measured amplitude 
is compared to its assumed functional form [Ref 6]: 
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( 1 ) 


y = A! Vr *e 

Dividing the amplitude at r by that at a range ro gives: 


y = y(ro)^^fr7Tr^e ( 2 ) 


Rearranging and taking the natural logarithm results in 


ln(V^ / ^ 0 * y I yo) = -oc(r - ro) (3) 

The slope related to the value of the attenuation coefficient, has a 
value of 2.44 dB/m which compares with results for similar shear 
waves derived by N. P. Chotiros, [Ref 8], for sand, Figure (3.3). By 



Figure (3.2) Value of attenuation coefficients. 
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plotting the data in this form, we can see that data should fail on a 
straight line with slope -a. The results for two different 
frequencies are shown in Figure (3.2). 



Figure (3.3) Attenuation measurements 
by Chotiros. 

B. AZIMUTHAL DEPENDENCE 

Having investigated the radial dependence of the wave 
amplitude, the next step is to investigate the azimuthal dependence. 

1. Procedure 

To insure relatively consistent conditions, the tank was 
watered and allowed to settle for 30 minutes to improve moisture 
uniformity. The source was put in place, and measurements were 
taken as rapidly as possible. Measurements were taken at one meter 
range with thirty degree seperation between locations. Data was 
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acquired using the vi illustrated in Figure (3.4), and a programing 
diagram in Appendix (B). This vi was designed to receive voltage 
amplitude triggered on the outgoing pulse of the HP3314A function 
generator and average a specified number of pulses to improve 
consistency in measurements. The results, for 400 Hz and 2000 Hz, 
shown in Figure (3.5), show that the spreading, was roughly uniform 
in all directions. It is interesting to note that these measurements 
were taken utilizing a one cycle pulse. Increasing the pulse length 
showed greater consistency in the received signal amplitude but 
created a larger problem with reflection off the walls of the tank. 
The results of the measurements, though not ideal, were considered 
adequate to attempt beamforming utilizing an array. 

C. SCATTERING FROM A BURIED OBJECT 

An evaluation of scattering from a buried object was 
conducted using a single source and multiple receiver set-up to 
evaluate the feasibility of using an array. This was accomplished by 
placing the source at the center of the tank and three receivers at 
approximately one meter ranges at various bearings out from the 
source. A vi. Figure (3.6), with a diagram and explanation in 
Appendix (C), was designed to acquire the received signal triggered 
on the source pulse and average it over a specified number of pulses. 
The averaged received signal was then saved, an a ten centimeter 
diameter, three centimeters thick, copper cylinder object placed in 
the tank covered by two centimeters of sand. The procedure was 
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repeated and the second averaged signal was subtracted from the 
first using the vi in Figure (3.7), with a diagram and explaination 
included as Appendix (C). As can be see in Figure(3.6), the arrival of 
the direct path and the reflected path are revealed in the time 
domain plot. By knowing the location of the source and the receiver, 
the speed of the Rayleigh wave can be determined from the direct 
path. The length of the relected path can then be determined. This 
length specifies an ellipse, corresponding to possible locations of 
the object. Figure (3.8). By utilizing multiple receivers, an object 
can be located by finding the intersection of multiple ellipses. This 
was accomplished with the data shown in Figures (3.9, 3.10, 3.11) 
and plotted in Figure (3.12). The location of the source and three 
receivers are indicated by the open circles. The location of the 
buried object is indicated by the x. It is seen that the x is near the 
interections of the ellipses. This techique must be developed 
further but the results are encouraging. 
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Figure (3.5) Azimuthal dependence of waveform 
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Figure (3.10) Comparison of 036.2 degree path. 
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IV. PHASED ARRAY BEAMFORMING 


Having investigated the single point surface wave source, 
attention is now turned to an array of sources. 


A. BEAMFORMING 

To investigate the feasibility of beamforming required some 
changes to the lab set-up. A three-element array configuration was 
selected, which required two additional sources to be built. Because 
the NB-MIO-16XL DAQ card can only support two analog output 
channels, a NB-AO-6 analog output card was acquired providing 
additional analog output channels. This addition enabled the use of a 
multichannel analog output control vi for the array. 


B. BEAMFORMING MODEL 

To predict the beam pattern, a model was derived 
based on the simple line array [Ref 5]. A line array of N simple 
sources spaced by a distance d, having the same amplitude and phase 
will generate a velocity of the form 


V{r, e,t) = 



j(wt-krt) 


( 4 ) 


where r/ is the distance from the /^^source to the receiver at 
{r,6). A program written in MATLAB, Appendix (D), produced a plot 
of the beam pattern based on the number of elements, their spacing. 
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amplitude weighting and phase. The results of this program are 
compared with the results acquired in the sand tank. 

C. EXPERIMENT 

The experiments were designed so as to allow comparison with 
the beamforming model. Phasing and weighting for each element 
were identical in the inital experiment. To evaluate the ability to 
beamform, the three elements of the array were fixed on the cross 
beam while the receiver was free to move azimuthal at fixed range 
from the center element. The Beamforming vi utilized, Figure (4.1), 
diagram Appendix (E), allowed for independent control of each 
element in frequency, phase, amplitude and pulse width to a high 
degree of accuracy. Though not utilized to its full capacity during 
these intitial experiments, it is able to support more complex array 
analysis. 

In the first experiment, the spacing was set to one 
wavelength {X = .18 meter) having measured a corresponding 
surface wave speed (c = 61.2 m/s) at a frequency of 340 Hz. The 
drive voltage and phase were equal for all elements. The pulse 
length was adjusted to fifteen cycles to allow for sufficient overlap 
of the signals from the three sources to produce interference. 

Figure (4.2) shows the received waveforms at eight selected 
bearings, arranged in approximately the correct orientation. The 
center plot represents a typical waveform received at bearings 
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3-Element Radiation Pattern for Theoreticai/Experimental data 



270 

340 Hz, 15 cycle pulse 


Figure (4.3) Calculated beam pattern for three element 
array overlayed by experimental data. 
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centered between the other locations. The variation in amplitude 
from bearing to bearing, and the symmetry about the axis 
demonstrates the beamforming. When compared to the theoretical 
beam pattern, Figure (4.3) it is seen that the array is producing the 
correct radiation pattern. Figure (4.4) shows the effect of 
beamforming at a specific bearing. The upper figure is the signal 
from a single source, the lower is the signal with all three sources 
driven. Clearly beamforming works. 
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V. CONCLUSION 


A. RESULTS 

The results of this research are encouraging. We successfully 
demonstrated target locatization using surface wave scattering and 
beamforming with phased array technology. Further development of 
the system is recommended. 

The size of the test tank proved to be a significant limitation. 
A next step would be to build a portable system that could be taken 
to the beach. Performing measurements at the beach will allow the 
use of longer ranges and lower frequencies. It will also provide 
realistic background noise. 

Other recommendations are to investigate pulse shaping and 
more advanced signal processing, techniques to address the issue of 
probability of detection, etc. This project is an ideal area to apply 
recent advances involving time reversal signal processing, [Ref 6]. 
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APPENDIX A. FACTSHEET ON SOURCE ELEMENT 


A four inch woofer speaker manufactured by Radio Shack was selected 
as a source element. Keeping cost considerations in mind, it had 
excellent source element characteristics. The long cone excursion of 
flexible rubber provided the necessary pressure to maintain coupling 
with the ground. The dust cap was sturdy enough to support various 
attachments affixed without deforming. The backplate was a large 
flat area to allow for attachment of a mounting post, Figure (A.1). 


WOOFER SPEAKER 
4'(10cm) woofer 

FEATURES 

• Smooth trequeocy response 

• Light weight hard cone paper based on careful 
study of internal toss values 

• Uproiled rubber edge minimizes distortion and 
improves linearity 

• Long travel'voice coil for high power handling- 



SPECIFICATIONS 

Impedance : e ohms 

Frequency response : 50-7.000Hz 

Free air resonance : SSHz+SHz 

Ots . 0.35 

Vas 0.23 eft 

vas is the volume of air that Is acoustically 
eqaivalent to the compliance of the cone 
SPL ; 84+2dB/W (1m) 

Input power nominal) ; SW 

Input power (maximum) : tOW 

Magnet weight : 228 g 8.0 oz 

Speaker : 700 g 24.7 oz 



i: 

li 

5 

li 

;! 

il 


il 


HOW TO USE 

1. Installation: 

a install the woofer unit from the enclosure front and keep it flush with the baffle surface, 
b Secure the speaker so that no air leakage is permitted. 

2. Connections.when using with a tweeter in a two-way system, connect as shown; 



Note When using the 40-1296 Crossover choose the 4000 Hz crossover point. 


Figure (A.1) Factsheet on source element 
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APPENDIX B. LABVIEW VI FOR INTENSITY FIELD STUDY 


The purpose of this vi is to take field intensity measurements 
to determine attenuation and to measure spreading. The vi is 
designed to be triggered by the HP 3314A Function Generator, 
acquire a specific number of received wave forms and average them. 
The three programing elements necessary were a "for loop," shift 
registers and an intermediate vi call Al wave. The three 



Graph or fiiel 


Figure (B.1). Elements of Intensity Field vi. 


elments can be seen in the simplified version above, Figure (B.1) and 
as elements of the actual program in Figure (B.2). The book is a "for 
loop." The arrows on each side are the shift registers. The Al WAVE 
vi which acquires the signal is also shown. In this case the "for 
loop" value is 100. Therefore, 100 aquisitions are made, summed 
and then divided by the "for loop" value to get an average of the 
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signal. The shift registers are used to carry the summed value to 
the beginning of the loop until the "for loop" value is reached. This 
is a timed acquisition, meaning that a hardware clock is used to 
control the acquisition rate. It is also a buffered acquisition, 
meaning that the data is stored in an intermediate memory buffer 
after it is acquired from the analog input channel. On the left side 
of the vi, in the lettered boxes. Figure (B.2), are the front panel 
inputs controlling the number of iterations, device, channel input, 
trigger type, and scan information. 
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Figure (B.2) Field intensity vi. 
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APPENDIX C. LABVIEW VI'S FOR DATA ACQUISITION 


Two vis were used to acquire and analyze data to determine 
object location in the sand by surface wave sattering. The initial 
aquisition vi was similar to the vi used to take field intensity 
measurements except it created a binary file that could be used by a 
comparision vi. To create this file an addition was made to the field 
intensity measurements vi as shown below, Figure {C.1). This can be 
seen in the vi on the far right of Figure (C.2). 

The comparision vi was written to read two binary files, 
displaying the first in the upper graph, both sets of points in the 
second graph, and the difference of the two in the third graph. In 
this way the additional disturbance in the field could easily be 
located in the time domain as shown in earlier Figures 
(3.9,3.10,3.11). 

Once the arrival times for the direct and relected paths had 
been determined, a plotting program written in MATLAB was created 
to plot the ellipses with a known distribution. This allowed for 
generating points having the highest probability for object 
detection. 
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Figure {C.1). Elements of Data Acquisition vi. 



Figure (C.2) Data Acquisition vi. 
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APPENDIX D. MATLAB PROGRAM FOR ELLIPSE 

PLOTTING 

As shown in Figure (3.8), it is possible to create a set of 
points that represent all the possible location of the target given 
the arrival time of the reflected and directed path. Using MATLAB, a 
program was developed to generate a set of ellipses representing 
multiple receivers. The program relies on the inputs provided by 
the user shown seen in the first few lines below: 

%Elipse Plotting of Scattering 

%Create it so the input radial angle is the axis for the elipse 

%Angle can be input converted to x, y coordinates 

home 

clear 

clg 

%Receivers 

for R=1:3; %total number of receivers 

bearing = input ('Enter the beamg to the Receiver ’); 
range = input ('Enter the range to the Receiver ’); 

TT = input ('Enter Travel Time of Reflected Path '); 

T = input ('Enter Travel Time of Direct Path '); 

A = range*sin((pi/180)*bearing); 

B = range*cos((pi/180)*bearing); 

xo = A/2; %center of elipse 

yo = B/2; %center of elipse 

delT = TT-T; %difference in travel times 

foe = range/2; %center to foci distance 

C = range/T; %wave speed based on direct path 

D = C*TT; %distance traveled by reflected path 

by = sqrt {(DI2)'^2 - foc'''2); %minor axis 

ax = D/2; %major axis 

t=0:.01:2*pi; 
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a=ax*cos(t); %ax is 1/2 length of elipse 

b=by*sin(t); %by is l/2width of elipse 

X=(a)+xo; %xo is center in x dir 

Y=b; %yo is center in y dir 

Z = X+Y*i; %z=X+Y*i 

Zp(R,:) = Z * exp(i*(pi*(0.5-(bearing/180)))); 


end 

clc 

dispCPick Mine Coordinates off the graph!') 

pause(5) 

figure(l) 

hold on 

for R = 1:3; % Change for the # of receivers 

plot(Zp(R,:)) 
end 

plot(0,0,'o') 

plot(76,-6,'o’) 

plot(77,10.3,’o') 

plot(52,-38,’o’) 

plot(86,-3i;x') 

hold off 

text(2,2,’source') 
text(80,-10, ’#1’) 
text(70,35, '#2') 
text(50,-50, ’#3’) 
text(47,-22,'object') 
xlabel(’range (cm)’) 
ylabel('range (cm)') 
grid 

[K,L] = ginput(l) 
clc 

disp(['Mine is located at X = ',num2str(K),' and Y = ',num2str(L)]) 
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APPENDIX E. MATLAB PROGRAM FOR BEAM FORMING 

ARRAY 


To accurately evaluate the three element array in it's ability 
to successfully beamform, a MATLAB program was written, with the 
help of Professor A. Atchley, Naval Postgraduate School. Given 
inputs of the distance to the receiver, spacing between elements, 
wavelength, and phase and amplitude for each element it accurately 
predicts the radiation pattern of the array in the XY plane. The 
following is the actual program listing: 

% Plot program for 3 element beam forming in the far field 


theta = (-pi:.01:pi); 

%angles 

r=l; 

%Distance to receiver from ce 

d = .18; 

%Spacing between elements 

kd = 0.9.*2*pi; 

%wave number * d 

k = kd/d; 

%wave number 

% Source 1 is at the top, 3 is at the bottom. 

fl = 0*pi/3; 

%phase control of element 1 

f2 = 0*pi/3; 

%phase control of element 2 

f3 = -0*pi/3; 

%phase control of element 3 

al = 1; 

%amplitude of element 1 

a2 = l; 

%amplitude of element 2 

a3 = l; 

%amplitude of element 3 

% Range from element to receiver 

rl = sqrt(r*r*cos(theta).*cos(theta) + (r*sin(theta) - 

d).*(r*sin(theta) 

-d)); 

r2 = r; 
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r3 = sqrt(r*r*cos(theta).*cos(theta) + (r*sin(theta) + 
d).*(r*sin(theta) + d)); 

% Pressure created by each element 

pi = al *exp(i*(-k*rl + fl))./sqrt(rl); 
p2 = a2*exp(i*(-k*r2 + f2))./sqrt(r2); 
p3 = a3*exp(i*(-k*r3 + f3))./sqrt(r3); 

pabs = abs(pl+p2+p3); 

pnorm = pabs./niax(pabs); 

%Experimental Data 

Rdata = [380,110,110,110,130,130,150,150,380,380,280, 
150,100,... 

150,150,50, 80,180,360,280,100,100,120,100, 80, 80, 

80, 340,... 

280,140, no, 130,160,130,140,150, 380]; 

ttheta=0 :pi/l 8:2*pi; 
data = Rdata/380; 
polar(theta,pnorm); 
hold on 

xlabel(’375 Hz, 15 cycle pulse’) 

title(’3-Element Radiation Pattern for Theoretical/Experimental 

data'); 

grid; 

polar(ttheta,data) 
hold off 
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APPENDIX F. LABVIEW VI FOR BEAMFORMING ARRAY 


A LABVIEW vi was written to control the three element array. 
The basic design involved creating a sine wave of one cycle with a 
specific amplitude and frequency inside a "for loop" which generated 
multiple cycles. This was then multiplied by a square wave that had 
been created in a similar fashion. The results is a pulsed sine wave. 
These elements can be seen in the simplified version, Figure (F.1). 



Figure (F.1). Elements of Three Element Function Generator vi. 


This basic process is repeated three times, one for each element 
and then combined into one matix with each row representing the 
output of one function generator. This matrix is then sent to an 
analog output write file which repeats the data inside a "while loop" 
until it is told to stop. Figure (F.2) is the diagram of the working 
program. 
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